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A novel, efficient air data calibration method is proposed for aircraft with 
limited envelopes. This method uses output-error optimization on three-dimensional 
inertial velocities to estimate calibration and wind parameters. Calibration 
parameters are based on assumed calibration models for static pressure, angle of 
attack, and flank angle. Estimated wind parameters are the north, east, and down 
components. The only assumptions needed for this method are that the inertial 
velocities and Euler angles are accurate, the calibration models are correct, and that 
the steady-state component of wind is constant throughout the maneuver. A 
two-minute maneuver was designed to excite the aircraft over the range of air data 
calibration parameters and de-correlate the angle-of-attack bias from the vertical 
component of wind. Simulation of the X-48B (The Boeing Company, Chicago, 
Illinois) aircraft was used to validate the method, ultimately using data derived from 
wind-tunnel testing to simulate the un-calibrated air data measurements. Results 
from the simulation were accurate and robust to turbulence levels comparable to 
those observed in flight. Future experiments are planned to evaluate the proposed 
air data calibration in a flight environment. 
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Nomenclature 

standard sea level speed of sound 
differential pressure bias 
differential pressure inverse gain 
differential pressure flank angle gain 
flank-angle influence on angle-of-attack gain 
angle-of-attack influence on flank angle gain 
upwash 
sidewash 

total temperature recovery factor 
Mach number 
total pressure 

calibrated static pressure 
measured static pressure 
calibrated differential pressure 
measured differential pressure 
ambient temperature 
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total temperature 

standard sea level temperature 

airspeed component in the body fixed x direction 

true airspeed 

computed groundspeed in the east direction 

computed groundspeed in the down direction 

computed groundspeed in the north direction 

windspeed in the east direction 

windspeed in the down direction 

windspeed in the north direction 

airspeed component in the body fixed y direction 

airspeed component in the body fixed z direction 

flank-angle bias 

calibrated flank angle 

measured flank angle 

angle-of-attack bias 

calibrated angle of attack 

measured angle of attack 

calibrated angle of sideslip 

ratio of specific heats for air 

pitch angle 

bank angle 

heading angle 


I. Introduction 

S imultaneous calibration of static pressure, angle of attack, and flank angle over the range of the flight 
envelope is possible with appropriate maneuver design and analysis methods. Air data accuracy is 
critical for pilot displays, as feedback to flight control laws, and for use in post-flight data analysis. A 
time-efficient method, appropriate for vehicles with a limited envelope, is presented and tested in a 
nonlinear simulation. This method uses a dynamic maneuver to excite the aircraft over a range of airspeed, 
angle of attack, and sideslip. Analysis is performed using output-error optimization and assumed calibration 
models. 

Static air data calibration methods involve maneuvers at many indicated airspeeds, often with an 
outside observer. Some of these methods include: control tower flybys, trailing cone with pacer aircraft, 
and radar tracking. 1,2 Some additional static methods that take advantage of a global positioning system 
(GPS) include three-leg, 2 3 four-leg, and cloverleaf maneuvers, 4 and air data circles. 5 Dynamic maneuvers 
have been proposed and used in flight to calibrate static pressure over a range of airspeeds. 6,7 The proposed 
method in this paper makes extensions to this dynamic maneuver to include calibrating angle of attack and 
flank angle while taking into account vertical wind. 

This proposed air data calibration method was tested in the X-48B (The Boeing Company, Chicago, 
Illinois) nonlinear simulation and included gradually increasing fidelity. The X-48B envelope is limited to 
incompressible airspeeds and altitudes below 10,000 ft above mean sea level. Initially, pre-determined 
calibration errors were introduced without turbulence for two wind conditions to check the accuracy of the 
estimated parameters and wind speeds against the true values. Then the same cases were run with 
turbulence levels representative of those observed in flight to test the robustness of the air data calibration 
method. Finally, simulated un-calibrated measurements, based on wind-tunnel data, were used to present 
the air data calibration algorithm with data that did not exactly match the assumed calibration models. This 
paper describes the aircraft and air data system, the proposed air data calibration method, and the results of 
the simulation testing. 
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II. Aircraft Description 

The X-48B aircraft, shown in Fig. 1, is a joint partnership between the National Aeronautics and Space 
Administration (NASA), the Air Force Research Laboratory (AFRL), and The Boeing Company. Built by 
Cranfield Aerospace Limited (Bedford, United Kingdom), the aircraft is a remotely-piloted, 8.5-percent 
dynamically-scaled hybrid-wing-body (HWB) aircraft with two air data probes; each probe includes 
pitot-static sensors, one angle-of-attack vane, and one flank-angle vane. 



Figure 1. The X-48B aircraft. 

The angle-of-attack and flank-angle vanes use optical encoders; both vanes were initially calibrated in 
the wind tunnel to account for local flow around the vehicle. The locations of the air data probes are such 
that the static pressure and the angle-of-attack vanes are affected by sideslip, and the flank-angle vanes are 
affected by angle of attack. Data tables from the wind-tunnel calibration are used within the X-48B 
nonlinear simulation to simulate un-calibrated static pressure, angle of attack, and flank angle. 

During the X-48B flight research program, eight three-leg maneuvers were performed early in the 
flight program to calibrate the air data. Pitot-static data were further calibrated with seven wind circles. 


III. Air Data Calibration 

The following three subsections describe output-error parameter estimation applied to air data 
calibration, the equations used for the calibration and output models, and the maneuver performed in 
simulation. 

A. Output-Error Parameter Estimation 

Output-error parameter estimation is the basis for the proposed air data calibration method. A diagram 
of the output-error parameter estimation processes is contained in Fig. 2 (Ref. 8). Inputs in this case are: 
total pressure, static pressure, total temperature, angle of attack, flank angle, pitch angle, bank angle, and 
heading angle. Outputs are the north, east, and down components of groundspeed. 
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Figure 2. Output-error parameter estimation (Ref. 8). 

Models are created for the calibration of static pressure, angle of attack, and flank angle based on a set 
of parameterized biases and gains. Initial guesses are supplied to the parameters. Parameters are updated 
based on the sensitivity of the model output to a change in the parameter. These parameters are used in a 
model to compute the north, east, and down components of groundspeed. Iterations are performed while 
updating the parameter estimates to minimize the residual. The residual is the error between the measured 
components of groundspeed and the estimated components of groundspeed. If the calibration model 
accurately describes the physical aircraft, the final residuals should be white noise. Reference 8 contains 
more information on parameter estimation and how to determine the model structure. 


B. Model Equations 

Equations are used to calibrate the air data and generate a computed groundspeed in the north, east, 
and down directions. These equations require the estimated parameters, presented in Table 1, and the 
sensor measurements, presented in Table 2. 

Table 1. Estimated parameters. 


Estimated parameters 

Description 

K\ 

Differential pressure bias 

K 2 

Differential pressure inverse gain 

K 3 

Differential pressure flank-angle gain 

K 4 

Flank-angle influence on angle-of-attack gain 

^5 

Angle-of-attack influence on flank-angle gain 

K a 

Upwash 

*a F 

Sidewash 

a b 

Angle-of-attack bias 

a Fb 

Flank-angle bias 

V 

r wn 

Windspeed in the north direction 

v we 

Windspeed in the east direction 

v wd 

Windspeed in the down direction 
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Table 2. Sensor measurements. 


Sensor measurements 

Description 

p T 

Total pressure 

P; 

Static pressure 

T t 

Total temperature 

a z 

Angle of attack 

°Fz 

Flank angle 

0 

Pitch angle 

O 

Bank angle 

W 

Heading angle 


Un-calibrated differential pressure is defined by the difference between total pressure and un-calibrated 
static pressure, Eq. (1). Differential pressure is calibrated using Eq. (2). The calibration is modeled as a bias 
plus a gain based on the inverse of un-calibrated differential pressure as recommended in Ref. 6 and 
depicted in Fig. 3. During simulation testing, it was noticed during analysis that the K 2 parameter should 
not be included in the calibration model. The parameter was removed from the model and not analyzed for 
this paper, but is included in Fig. 3 and Eq. (2) for the sake of completeness and general applicability. This 
equation is suitable for incompressible flow and an additional parameter may be required if compressibility 
is of concern. For the X-48B simulation, a gain on flank angle was added due to the effect of sideslip on 
static pressure. 



Figure 3. Example airspeed calibration curve shape. 

Total pressure is assumed to be correct (Ref. 2), so the static pressure is calibrated with the calibrated 
differential pressure, Eq. (3). 


A P Z =P T -P Z 


( 1 ) 
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P c = P T - AP C (3) 

Calibrated static pressure is used to compute Mach number, Eq. (4). Mach number and total 
temperature are used to compute temperature, Eq. (5), and true airspeed, Eq. (6). Due to the limited 
airspeed and altitude envelope of the X-48B aircraft, the total temperature recovery factor, k , was assumed 
to be unity and all of the error in airspeed was collected in the static pressure term. Level 
acceleration/deceleration maneuvers could be performed to independently estimate the total temperature 
recovery factor (Ref. 2). Using the approach of Eqs. (1) through (6), the total temperature recovery factor 
could not be treated as an additional estimated parameter because of the strong correlation, 0.96, between 
static pressure and temperature. 
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V = 



( 6 ) 


Angle of attack and flank angle were calibrated using Eqs. (7) and (8), which were modeled with bias 
and gain terms based on Ref. 9. Due to the locations of the vanes on the X-48B aircraft, a gain on flank 
angle was added to the angle-of-attack calibration and a gain on angle of attack was added to the 
flank-angle calibration. 


n , a z- a b xVfr 

°c - — + K 4 a Fz 
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K. 


L 5 u z 
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(7) 

( 8 ) 


Flank angle was converted to angle of sideslip using Eq. (9) (Ref. 2) and the lateral (Eq. 10), forward 
(Eq. 11), and down (Eq. 12) velocities were computed. 


tan 1 (tan(« Fc )cos(« c )) 

(9) 

v = V sin (/3 C ) 

(10) 

= V cos(a c )cos(/? c ) 

(11) 


6 

American Institute of Aeronautics and Astronautics 



w = V sin(or c )cos(/? c ) 


( 12 ) 


The aircraft velocity components and inertial angles are used with estimates of the forward, lateral, and 
down windspeed to compute north, east, and down inertial velocity, Eq. (13). These computed inertial 
velocities are compared with measured inertial velocities as part of the output-error parameter estimation. 


'cos(w)cos(0) (cos(w)sin(0)sin(o)-sin(w)cos(<X>)) (cos(w)sin(0)cos(o)+ sin(w)sin(<I>))' 

V 


v wn 

sin(w)cos(©) (sin(w)sin(0)sin(o)+ cos(w)cos(<X>)) (sin(w)sin(©)cos(o)- cos^^in^)) 

V 

+ 

V we 

-sin(0) cos(©)sin(<I>) cos(©)cos(<X>) 

w 


v wd 


(13) 


Finally, in addition to the model structures for calibration, it was assumed that the inertial angles and 
velocities are accurate and the steady-state wind remains constant throughout the maneuver. 

C. Maneuver Design 

Three goals of the maneuver design are to: estimate the wind, fully explore the ranges of estimated 
parameters, and de-correlate the estimates. It was desired to turn the aircraft through 180 deg of heading 
change in order to estimate the north and east winds. Airspeed, angle of attack, and flank angle must be 
excited in order to estimate their calibration parameters. Angle-of-attack bias and the vertical wind 
parameter are very nearly correlated; to de-correlate these estimates, bank angles near 45 deg were used 
during a portion of the maneuver. 

The maneuver design used for this simulation study is shown in Figs. 4(a) through 4(d). Due to 
proprietary considerations, the figures have been normalized. The maneuver consists of a shallow bank turn 
through 180 deg of heading change. During this turn, a power level angle (PLA) ramp reduces the airspeed, 
Fig. 4(a), while altitude is maintained, resulting in an increase in angle of attack, Fig. 4(b). This portion of 
the maneuver results in an estimate of north and east wind, static pressure calibration over the range 
excited, and angle-of-attack calibration over the range excited. A shallow bank angle was chosen to 
maintain quasi-steady flight conditions. Separating the turn from the airspeed reduction would be warranted 
if there were a single air data probe located off centerline rather than averaging dual air data probes or 
taking into account the yaw rate. Feft and right air data probes were averaged in this simulation study. 



Figure 4(a). True airspeed. 
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angle of attack 







Figure 4(d). Angle of sideslip. 

Following the turn, the PL A is returned to its cruising level and the aircraft is banked approximately 
45 deg and back to wings-level, Fig. 4(c). The bank is added to reduce correlation between the vertical 
component of wind and angle of attack. Table 3 is a comparison of correlation without and with the bank 
component of the maneuver; values closer to 1 indicate higher correlation. 

Table 3. Effect of 45-deg bank on correlation of angle of attack and vertical wind. 



Without bank 

With bank 

Correlation 

0.75 

0.57 


This maneuver is followed by a rudder pedal step, Fig. 4(d), in order to estimate the flank- angle 
calibration over the range of sideslip excited. The total time to complete the maneuver is approximately two 
minutes. 


IV. Results 

Air data calibration maneuvers were flown in the X-48B nonlinear simulation with two different wind 
cases, with and without turbulence; fidelity was gradually increased. A summary of the tests follows: 

1) Pre-determined calibration errors without turbulence 

2) Pre-determined calibration errors with turbulence 

3) Wind-tunnel calibration errors with turbulence. 

Wind Case 1 was a 13-kn wind out of 335 deg with an angle of 6 deg above the horizon. Case 2 was an 
8-kn wind out of 135 deg with an angle of 3 deg above the horizon. 

Pre-determined calibration errors were introduced to simulate un-calibrated data that exactly matched 
the model structure given in Eqs. (2), (7), and (8). Additionally, the X-48B nonlinear simulation includes 
air data calibration tables generated from wind-tunnel and flight data. These calibration tables allow the 
un-calibrated data to be used as input to the parameter estimation and provide the true value to validate the 
parameter estimation results. The turbulence levels used were representative of those observed in flight. 
The output-error optimization routine in the System IDentification Programs for AirCraft (SIDPAC) 9 10 
software was used, along with Eqs. (1) through (13), to estimate the calibration parameters and 
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three-dimensional wind. These parameters were used to calibrate the simulated time history data of static 
pressure, angle of attack, and flank angle. 


A. Estimation with Pre-Determined Calibration Errors 

Pre-determined errors were introduced in order to de-calibrate the equations of motion values from the 
simulation. These de-calibrated signals were then used as measurements in the air data calibration method 
to check its accuracy. Table 4 summarizes the results for wind Cases 1 and 2. Parameters K 3 , K 4 , and K 5 
were not included in the pre-determined calibration errors because they were based on each other, and an 
analytical solution to de-calibrate the static pressure, angle of attack, and flank angle with those parameters 
included could not be determined. These parameters were, however, included in the air data calibration 
with wind-tunnel-based data. 


Table 4. Case 1 and Case 2 results with no turbulence. 


Parameter 

Case 1 
true value 

Case 1 estimate 

Case 2 
true value 

Case 2 estimate 

K i 

7.00e-2 

7.09e-2 

7.00e-2 

7.08e-2 

K a 

1.60 

1.70 

1.60 

1.70 


1.05 

1.05 

1.05 

1.05 

a b 

1.20° 

1.66° 

1.20° 

1.66° 

a Fb 

0.60° 

0.61° 

0.60° 

0.61° 

V 

r wn 

-1 1.72 kn 

-11.69 kn 

5.65 kn 

5.67 kn 

v we 

5.46 kn 

5.47 kn 

-5.65 kn 

-5.65 kn 

V wd 

1.36 kn 

2.16 kn 

0.42 kn 

1 .22 kn 


From Table 4, it is apparent that the majority of parameters were estimated very accurately. 
Angle-of-attack bias and the vertical component of wind are the two parameters that were not estimated 
accurately, but accuracy was improved with the addition of turbulence, as discussed below. 

B. Performance of Air Data Calibration in the Presence of Turbulence 

Turbulence of a similar magnitude to that seen during flight-testing was introduced into the same two 
Cases discussed above. The results for Case 1 are shown in Table 5 and for Case 2 in Table 6. 


Table 5. Case 1 results with turbulence. 


Parameter 

Case 1 
true value 

Case 1, 
no turbulence 

Case 1, 

with turbulence 


7.00e-2 

7.09e-2 

6.87e-2 

K a 

1.60 

1.70 

1.71 

*a F 

1.05 

1.05 

1.06 

a b 

1.20° 

1.66° 

1.33° 

a Fb 

0.60° 

0.61° 

0.58° 

V 

r wn 

-1 1.72 kn 

-1 1.69 kn 

-1 1.74 kn 

v we 

5.46 kn 

5.47 kn 

5.46 kn 

v wd 

1.36 kn 

2.16 kn 

2.00 kn 
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Table 6. Case 2 results with turbulence. 


Parameter 

Case 2 
true value 

Case 2, 
no turbulence 

Case 2, 

with turbulence 

K i 

7.00e-2 

7.08e-2 

6.66e-2 

K a 

1.60 

1.70 

1.75 

«a r 

1.05 

1.05 

1.07 

a b 

1.20° 

1.66° 

0.92° 

a Fb 

0.60° 

0.61° 

0.55° 

V 

r wn 

5.65 kn 

5.67 kn 

5.58 kn 

v we 

-5.65 kn 

-5.65 kn 

-5.64 kn 

v wd 

0.42 kn 

1.22 kn 

0.96 kn 


Calibration parameter estimates were not greatly affected by the addition of turbulence except for 
angle-of-attack bias. Both angle-of-attack bias and the estimate of vertical wind were more accurate after 
the addition of turbulence. 

From this simulation it was concluded that the proposed air data calibration technique produced 
reasonable air data calibration and wind estimates. The testing up to this point involved pre-determined 
calibration errors, which exactly matched the models assumed for the calibration technique. The next step 
was to use simulated un-calibrated data from the X-48B nonlinear simulation, which are based on 
wind-tunnel testing, and compare the estimates to the true values. Assumed calibration models will no 
longer perfectly match the actual calibrations. 


C. Estimation with High-Fidelity Calibration Errors 

Simulated un-calibrated signals from wind-tunnel data were used as input to the air data calibration 
method. A comparison of the simulated and estimated wind is given in Table 7. 


Table 7. Wind estimation results. 



Case 1 

Case 1 

Case 2 

Case 2 

Parameter 

true value 

estimate 

true value 

estimate 

V 

r wn 

-1 1.72 kn 

-11.75 kn 

5.65 kn 

5.58 kn 

v we 

5.46 kn 

5.46 kn 

-5.65 kn 

-5.63 kn 

v wd 

1.36 kn 

1.37 kn 

0.42 kn 

0.35 kn 


Estimation of the speed and direction of the wind produced impressive results in both cases with very 
little error between the estimate and the truth data. Estimating the vertical component of wind was much 
more difficult for two reasons. First, the magnitude is very small, approximately 1.4 kn in Case 1 and 
0.4 kn in Case 2. These small magnitudes make the impact of the vertical wind on the vertical groundspeed 
small and, consequently, harder to estimate. Second, as can be seen in Table 3, angle-of-attack bias and 
vertical wind are somewhat correlated, complicating estimation of the vertical wind. 
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D. Comparison of Air Data Calibration with Equations of Motion Values 

Air data calibration parameters for wind Case 1 are given in Table 8. Also included with the estimates 
are Cramer-Rao bounds for the parameters. 


Table 8. Case 1 parameter estimates. 


Parameter 

True value 

Estimate 

Cramer-Rao bound 

K i 

_ 

7.55e-2 

1.54e-5 

K 3 

_ 

-6.24e-4 

2.17e-6 

K 4 

- 

6.26e-l 

9.29e-3 

^5 

_ 

-1.36e-l 

4.06e-3 

K a 

- 

1.77 

1.75e-3 

*a F 

_ 

1.05 

5.66e-4 

a b 

_ 

2.46° 

6.07e-2 

a Fb 

- 

-2.05° 

4.99e-4 

V 

wn 

-1 1.72 kn 

-11.75 kn 

1.04e-3 

v we 

5.46 kn 

5.46 kn 

9.63e-4 

V wd 

1.36 kn 

1.37 kn 

4.87e-3 


The calibration parameters were used to generate calibrated time histories for the maneuver and are 
compared to the true values in Figs. 5(a) through 5(c) for true airspeed, angle of attack, and sideslip angle 
for wind Case 1. Wind Case 2 produced similar results. True values are from the equations of motion and 
auxiliary equations in the simulation. Figs. 5(a) through 5(c) show that the calibrated parameters match the 
true values quite well, indicating that the assumed calibration models were correct and the parameter 
estimates were accurate. 

North, east, and down groundspeeds used by the output-error optimization are included in Figs. 6(a) 
through 6(c). 



Figure 5(a). Wind Case 1 true airspeed calibration. 
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Normalized angle of sideslip 



Figure 5(b). Wind Case 1 angle-of-attack calibration. 



Figure 5(c). Wind Case 1 angle-of-sideslip calibration. 
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Normalized east groundspeed Normalized north groundspeed 
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Figure 6(a). Wind Case 1 north groundspeed optimization. 
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Figure 6(b). Wind Case 1 east groundspeed optimization. 


14 

American Institute of Aeronautics and Astronautics 



Figure 6(c). Wind Case 1 down groundspeed optimization. 

Figs. 6(a) through 6(c) show that the computed output was optimized to match the measured output 
very well. This observation indicates that at least the minimum number of model parameters was estimated 
and that the output-error optimization routine was working properly. 

Testing this proposed air data calibration method included gradually increasing fidelity. This process 
started with pre-determined de-calibration parameters and no turbulence. Parameter estimates in this case 
were accurate with the exception of angle-of-attack bias and vertical wind. Adding realistic levels of 
turbulence to the simulation improved estimation of angle-of-attack bias and vertical wind while retaining 
the accuracy of the estimates of the other parameters. Finally, a simulation was run with wind-tunnel-based 
de-calibration parameters and turbulence. This last step increased the fidelity by using an assumed 
calibration model that did not perfectly match the de-calibration. Results were the most accurate in this case 
with good de-correlation of angle-of-attack bias and vertical wind. From this series of simulation tests, it 
was concluded that the proposed air data calibration method is feasible for use in a flight environment. This 
test series also presents an efficient technique that does not require extensive airspace or flight time and that 
should be appropriate for aircraft that are equipped with a good inertial navigation solution. 

V. Conclusions 

A novel air data calibration method has been proposed and evaluated in simulation. This method 
estimates calibration parameters for static pressure, angle of attack, and flank angle over the range excited 
in addition to three-dimensional winds using a single maneuver. Output-error optimization is used on the 
north, east, and down inertial speeds. The only assumptions necessary are that the inertial speeds and Euler 
angles are accurate, the assumed calibration models are correct, and that the steady-state wind is constant 
during the maneuver. Simulations were run with two wind cases and turbulence using both pre-determined 
and wind-tunnel-based de-calibration parameters. Results from the simulation evaluation were that the 
proposed method produced accurate estimates of the parameters and successfully de-correlated 
angle-of-attack bias and the vertical component of wind. Results were accurate for the highest fidelity 
simulation run, which used data from wind-tunnel testing for the de-calibration parameters and included 
turbulence. 

The proposed air data calibration method has the potential to significantly reduce the time and the 
airspace required for air data calibration of aircraft. It is possible to perform calibration of static pressure, 
angle of attack, and flank angle over the entire envelope with a single, two-minute maneuver. This research 
and simulation testing was performed on a vehicle with a limited envelope and incompressible flow. Future 
research could be to expand the air data calibration method to larger envelopes, compressible flow, and 
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include simultaneous calibration of the total temperature recovery factor. Additionally, the amount of wind 
variation allowable during a maneuver and its effect on air data calibration accuracy could be investigated. 
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